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^ Abstract 

, — J In the AdS / CFT framework meson thermalization in the presence of a con- 

stant external magnetic field in a strongly coupled gauge theory has been 
studied. In the gravitational description the thermalization of mesons corre- 
sponds to the horizon formation on the flavour D7-brane which is embedded 
in the AdS§ x S 5 background in the probe limit. The apparent horizon forms 
due to the time-dependent change in the baryon number chemical potential, 
the injection of baryons in the gauge theory. We will numerically show that 
the thermalization happens even faster in the presence of the magnetic field 
on the probe brane. We observe that this reduction in the thermalization 
time sustains up to a specific value of the magnetic field. 
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1 Introduction and Results 

Quark-gluon plasma created by colliding relativistic heavy ions such as Gold 
at Relativistic Heavy Ion Collider or Lead at Large Hadron Collider, is be- 
lieved to be strongly coupled [1,2]. In the initial stages right after the collision 
the quark-gluon plasma is far from equilibrium and hence it would be dif- 
ficult to study its properties. But after a very short time, (r < 1 fm/c), 
the hydrodynamics is applicable indicating that the local equilibrium has 
been partly reached [3,4]. An important question is how a strongly coupled 
system can rapidly relax from far from equilibrium regime to the hydrody- 
namic phase. This process is generally called rapid thermalization. Since 
the system is strongly coupled the usual perturbative methods seem inap- 
plicable to study such a process. Recent calculations show that AdS/CFT 
correspondence [5] which loosely speaking is a duality between a strongly 
coupled system and a classical gravity theory can be a suitable framework to 
describe rapid thermalization [6-8]. 

An interesting discussion has been made in [9] indicating that after the 
collision of two heavy ions a strong magnetic field is produced which is present 
only for a very short time. Since this magnetic field is nonzero in the initial 
stages of the collision it is reasonable to study the process of thermalization 
in its presence. Therefore one can use the gauge/gravity duality to explain 
how the magnetic field can effect the rapid thermalization in QGP. This is 
exactly the question that we are trying to address in this paper. 

According to the AdS/CFT correspondence the four- dimensional super- 
conformal Yang-Milles (SYM) theory with gauge group SU(N C ) is dual to the 
type IIB string theory on AdS^ x S 5 geometry which describes the near hori- 
zon geometry of a stack of N c extremal D3-branes [5]. This correspondence 
is more understood in the large N c and large t'Hooft coupling, A = gy M N c , 
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limit. In these limits the type lib supergravity on AdS^ x S 5 geometry is dual 
to the strongly coupled SYM theory and hence it provides a useful tool to 
study the strongly coupled regime of the SYM theory. This correspondence is 
extended in various ways. In particular the thermal SYM theory corresponds 
to the supergravity in the AdS-Schwarzschild black hole background. The 
Hawking temperature of the black hole is identified with the temperature of 
the thermal gauge theory [5]. 

In the context of AdS/CFT correspondence matter fields (or quarks) in 
the fundamental representation are introduced in the gauge theory by em- 
bedding flavour D7-branes in the probe limit [10]. By probe limit one means 
that the number of flavour D7-branes are much smaller than the number of 
D3-branes. In this system open strings stretched between D3-branes and D7- 
branes are realized as quarks in the SYM theory. Their mass is proportional 
to the asymptotic separation of the D3 and D7-branes in transverse direc- 
tions. Additionally the open strings with both endpoints on the D7-branes 
are considered as mesons. This system is a notable candidate to describe 
QCD-like theories [11]. 

The process of thermalization refers to an increase in the temperature of 
a system from an initial value which can be zero to a higher value in a non- 
adiabatic way which takes the system out of equilibrium. This is done by the 
injection of energy into the gauge theory or equivalently, in the context of 
AdS / CFT, by the collapse of matter in AdS§ x S 5 background. This injection 
of energy can be done by turning on a time-dependent source term in the 
gauge theory. It has been shown that the addition of time-dependent sources 
such as metric or scalar field on the boundary result in horizon formation in 
the bulk [6,7]. These sources are coupled to the energy-momentum tensor 
and scalar operators, respectively. Therefore one can simply say that the 
thermalization in the field theory is dual to the black hole formation in the 
bulk. 

So far we have been discussing the thermalization in the so-called gluon 
sector which corresponds to the horizon formation in the bulk. This can be 
generalized to the horizon formation on the probe branes which is usually 
referred to as thermalization in the meson sector [12-14]. In [12] the energy 
injection that leads to the thermalization in gauge theory is modeled by the 
injection of baryons which corresponds to throwing the open strings on the 
probe D7-brane from the boundary. The time-scale of the thermalization 
is identified with the time-scale of the apparent horizon formation observed 
on the boundary. In such a set-up we turn on a constant magnetic field on 
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the D7-brane and study the meson thermalization in the presence of this 
magnetic field. 

The main result of this paper is that the thermalization time correspond- 
ing to the meson melting decreases as one increases the external magnetic 
field on the probe D7-brane. This can be justified by an energy evalua- 
tion which will be explained in more details in the last section. It has been 
observed that when the external magnetic field is zero the thermalization 
time scales as ( a A 2 ) 5 where Ub is the baryon number density which is in- 
jected into the system and u shows the period of time that the injection 
happens [12] . In the presence of the magnetic field the difference between 
the thermalization times of zero and non-zero magnetic fields behaves as 

where for weak magnetic field, defined as {2na')B <C 1, k is 2 and as B is 
raised it deviates from 2. This model shows that the presence of the magnetic 
field helps the thermalization, or in more exact words meson dissociation, to 
happen even faster. We observe that this reduction in the thermalization 
time happens up to a specific value of the magnetic field which depends on 
the numeric values of the parameters of the theory. 

2 D7-Brane Embeddings in External Mag- 
netic Field 

In AdS/CFT one can add the fundamental matter into the gauge theory 
by embedding a probe D7-brane in an asymptotically AdS background [10]. 
The asymptotic behaviour of this brane teaches us about the fundamental 
matter (quark) mass and the condensation which indicates the spontaneous 
chiral symmetry breaking [15]. One way to break the chiral symmetry is to 
embed the D7-brane in the AdS$ x S 5 background with a nonzero constant 
Kalb-Ramond field 3 [16-18]. This has the effect of repelling the D7-brane 
from the origin. The nontrivial shape of the probe brane gives the conden- 
sation (c) which varies with respect to the quark mass, c with respect to m 
shows a spiral behaviour near the origin. This means that for a given value 

3 Note that a pure gauge Kalb-Ramond field which is a solution to type lib supergravity 
equations of motion can be replaced by the constant magnetic field B on the D7-brane. 
This field content does not deform the AdS§ x S 5 background [17]. 
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of the quark mass a number of solutions exist. Among them the physical so- 
lution is the most energetically favourable one. The presence of the constant 
Kalb-Ramond or equivalently constant magnetic field leads to nonzero quark 
condensate even at zero quark mass. In the following we will discuss these 
results in more details. 

The background of interest is AdS^ x S 5 geometry describing the near 
horizon geometry of N c D3-branes 



ds 2 = {^) 2 {-dt 2 + dx 2 ) + (-)W + P 2 dn 2 3 + da 2 + a 2 d V 2 ) 
R u 

9s C (4) = -^dx° A dx 1 A dx 2 A dx 3 , g s = e 0o °, R A = Awg s N c l; , 
R 



(2.1) 



where Cu\ and are four-form and Dilaton fields, respectively. In this co- 
ordinate u 2 = p 2 + a 2 and l s = \foi' is the string length scale. In order to 
introduce the fundamental matter we have to add a space-filling flavour D7- 
brane in the probe limit to this background. The configuration of the D3-D7 
branes is 

0123456789 
D3 x x x x (2.2) 
D7 xxxxxxxx 

In the probe limit the low energy effective action for a D7-brane in an arbi- 
trary background is described by Dirac-Born-Infeld (DBI) and Chern-Simons 
(CS) actions 

S = 5'dbi + Scs ' 
5'dbi = — A*7 J d 8 £ e -0 y/- det(g ab + B ab + 27ra'F ab ) , ^ ^ 

5cs = »r J P[ZC {n) e B ]e 2 ™' F , 

where induced metric g ab and induced Kalb-Ramond field B ab are given by 

g a b = G MNd a X M d b X N , ^ 
B a b = B MN d a X M d b X N . 

£ a are worldvolume coordinates and fi^ 1 = {2n) 7 l s s g s is the D7-brane tension. 
Gmn is the background metric introduced in (2.1) and in our background 
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B MN is zero. F ab is the field strength of the gauge field living on the D7- 
brane. In the CS part, C( n ) is a ra-form field and P[...] denotes the pull-back 
of the bulk space-time tensors to the D7-brane worldvolume. 

We work in the static gauge where the D7-brane is extended along t, x, p, Q3 
and consider the magnetic field 

B = F xy , (2.5) 

where B is constant. The presence of this external magnetic field will break 
supersymmetry on the probe brane. a as a function of p describes the shape 
of the brane. We set all other fields to zero. The D7-barne action then 
becomes 

S DBI = -p 7 V 3 Vol(fi 3 ) J dtdp pV 1 + (2na') 2 R 4 u- 4 B 2 Vl + a' 2 , (2.6) 
where a' — The equation of motion for a(p) can be easily found 



/ pV / (2Tra') 2 B 2 R 4 \ 2p 3 (2W) 2 ^ 2 ^ 4 ayTTg 
"VVT+^V + (P 2 + ^ 2 ) 2 ) + (P 2 + ^ 2 ) 3 ,/i i (W)W ~ " 

(2.7) 

In the asymptotic limit where p — > 00, the above equation reduces to 

9,, 



pV 



and hence its solution approaches 

a{p) =m + 4 + -- (2.9) 
p 2 

According to the AdS/CFT correspondence the leading and sub- leading 
terms in the asymptotic expansion of various bulk fields define the source 
for the dual operator and the expectation value of that operator, respec- 
tively. For cr(p) the leading term corresponds to the quark mass, m q = 
m/(2ira'), and the sub-leading term is related to the chiral condensate via 

c= (27ra') 3 (qq) [15]. 

The equation of motion for <r(p) is complicated to be solved analytically 
and hence we use numerical methods to find its solutions. In order to solve 
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Figure 1: Left: The condensation as a function of the mass for 2na'B = 0.09 
and R — 1. Points M, N and P represent three massless solutions. Right: 
The shape of three massless solutions corresponding to the points M, N and 
P. The blue curve (representing point M) is the physical massless solution. 



it numerically we have to impose two boundary conditions at a specific point 
which we choose them to be <r(0) > and cr'(0) = 0. The second condition 
guarantees the regularity of the D7-brane embedding. The above boundary 
conditions for a fixed value of the magnetic field and different values of cr(0) 
leads to different values of mass and condensation. The final results are 
plotted in figure l(left). Evaluating the energy for different configurations 
of c and m shows that the region with m > and — c < 0, the bottom 
right quadrant in the figure l(left), describes physical solutions. It can be 
seen that in the presence of a fixed magnetic field there are three massless 
solutions. One of them is trivial and is zero for all values of p indicating that 
m = c = 0. Although for the other two solutions m is still zero, the value of c 
is not zero. A nonzero c means that chiral symmetry has been spontaneously 
broken by the external magnetic field [16]. Among massless solutions in 
figure 1 (left) the solution corresponding to the point M is energetically more 
favourable. cr(0) at this point is 0.139 where 2ita'B = 0.09. The blue, red and 
green curves in figure 1 (right) represent brane embeddings with asymptotic 
behaviours corresponding to the points M, N and P, respectively. 

In figure 2 physical massless configurations for various values of the mag- 
netic field are plotted. It is clear that for larger values of the magnetic field 
er(0) also becomes larger. In other words the turning point of the D7-brane 
configurations pushes more and more towards the boundary as one raises B. 

An important observation is that for small values of the magnetic field 
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the deformation of the D7-brane is not significant. Therefore a(p) ~ is an 
acceptable approximation in the weak magnetic field limit. We will use this 
fact in the subsection 3.1 to simplify the calculation of the thermalization 
time. 

dp) 

0.20 - 

0.15 

0.10 

0.05 

0.00 
-0.05 
-0.10 - 

Figure 2: Different physical massless solutions for (27ra')B = 
0.09, 0.07, 0.05, 0.03 and 0.01 (top to bottom) and R=l. 

3 Meson Thermalization in External Magnetic 
Field 

The process in which the temperature of a system changes non-adiabatically 
from Ti (which can be zero or non-zero) to a larger value T 2 is generally 
called thermalization. In the context of AdS / CFT a zero temperature state 
in a strongly coupled gauge theory is dual to the pure AdS$ x S 5 background 
in the bulk. Similarly a thermal state is dual to the AdS-bla,ck hole back- 
ground. Then the thermalization in the gauge theory side is identified by 
the horizon formation in the gravity side. A generalization has been done 
in [12] where meson thermalization is identified by the apparent horizon for- 
mation on the flavour D7-brane. The thermalization happens due to the 
injection of baryons into the gauge theory. In the gravity side baryon injec- 
tion is presented by throwing the fundamental strings on the D7-brane. The 
open strings stretched between D3- and D7-branes play the role of the fun- 
damental quarks in the gauge theory side. A time dependent non-adiabatic 
baryon injection will suddenly increase the number of baryon charges in a 
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time dependent manner and consequently produces a time dependent chem- 
ical potential. The time dependence of the chemical potential or baryon 
number density can finally lead to the horizon formation on the D7-brane. 

In the AdS/CFT correspondence the static chemical potential, p, is in- 
troduced via time component of the gauge filed as 

/■oo 

\x = \ dp dpAt(p), (3.1) 
Jo 

in the gauge where the radial component of the gauge field is set to zero. In 
order to obtain a time dependent chemical potential we have to consider a 
time dependent gauge filed configuration. This can be done by introducing 
source terms to the DBI action 

5S = /i 7 V 3 Vol(0 3 ) J dtdp(A t f + A p f), (3.2) 

which describes the dynamics of open strings, or equivalently, fundamental 
quarks injected on the brane from the boundary. Notice that the source 
terms are now time dependent and their form determines how the chemical 
potential varies in time. 

Since we consider the massless quarks it is convenient to work with light- 
cone coordinates which are defined as 



x ± 



t±z = tT J dpy/H(l + a'), (3.3) 



where H = (i?/w) 4 . Without loss of generality we assume the currents are 

only functions of x~ . More precisely we consider j p = —\H(1 + a')) 2 f = 
g'{x~). Therefore in the presence of the magnetic field equations of motion 
for A t and A p are 



(2W)F tp = 9VY +^ (3.4) 

P vV + (27ra') 2 p 6 (l + {2na') 2 HB 2 ) 

Notice that the above assumption for the currents satisfies current conserva- 
tion equation [14]. Moreover the charge associated with this current is 

Q= I d 7 ZV=g f, (3.5) 
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where Q is the number of quarks or equivalently open string endpoints on 
the D7-brane. Hence the baryon number density can be defined as 

hb = v^k' (3 - 6) 

and using (3.5) we finally have [14] 

g(x-) = (2/n)(2ira') 4 n B (x-). (3.7) 

Mesons in the AdS/CFT duality are realized as small fluctuations of the 
shape of the flavour D7-brane in the transverse directions [19]. In order 
to find the action describing the dynamics of scalar mesons we assume the 
transverse directions have the form 

x I = x I + y I , (3.8) 

where Xq are the classical solutions denoting the shape of the D7-brane and 
y 1 are the small fluctuations. In the presence of the magnetic field, B, and 
the gauge field strength, F tp , the action for mesons is obtained by expanding 
the DBI action for small fluctuations. Up to second order in y 1 the action 
becomes 

Sdbi = Sq + Si + ... 

= S - ^ [ d 8 ^ (< b M ba + ^ b N ba 



l -^ b M bc ^ Mda +I( 7o ^M ba ) 2 + ... 



(3.9) 



where 

70 ab 



G ab + Gudaxfaxl + (2W)F afc , (3.10a) 

M ab = djGaty 1 + Gud a xld h y J + Gjjd^dyxl + d K Gud a XQd h x J Q y K ', 

(3.10b) 

N ab = ) 2 didjGa b y I y J + Guday^by 1 + d K G IJ d a x I d b y J y K 

+ d K G I jd a y I d b xiy K + ^<9 L G 7J d a 4«9 b x J y V, (3.10c) 
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and 

S = -n 7 I d»^ . (3.11) 



One can easily observe that if B is set to zero in this set-up which means 
that the classical solution Xq is constant and not coordinate-dependent, the 
action in [12, 14] is recovered. It should be mentioned that the first order 
terms in y 1 are cancelled due to their equation of motion. Note also that the 
CS action does not contribute since the gauge field and the four-form field 
have common indices. 

In our case we consider the profile of the transverse directions as 

<T = *o(j>)+x(Z), (3-12) 

and ip = 0. We choose <7 (p) to be a massless solution of the equation of 
motion (2.7) and x is a small fluctuation around this solution. After some 
long calculations the part of the action which describes the kinetic term for 
the fluctuations in \ reduces to 



Si = ~ f J d 8 £ G aa [(1 + (W) 2 1T2J 2 ) (1 - (W) 2 f* + <x 2 ) 

x a/- det s §Q b d a xd b x, 

where we have defined 



(3.13) 



1 - {I-kcy'YF 1 

ab _ I " > tp ab a a Pa b Prr^r (1 \£\ 

- i /o A9Z72 , _/2 S _ ^0 y a ^va- (A i4 J 



1 - (2na>yF t l + a' 2 °° 



Sq 6 and 9q are the symmetric and anti-symmetric parts of 7q 6 , respectively 

7g 6 = «o 6 + C (3-15) 
In our set-up the only nonzero component for 0Q b is 



Hence the action S\ for the scalar mesons can be written more neatly as 
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6 

) 



^1 
(i 

) 



where the explicit form of the metric components become 

-S« = f4 (1 + (2na') 2 HB 2 )^ (l - (2W) 2 F 2 p + of) 

s pp = /if (1 + (2ira') 2 HB 2 )^ (l - (2rra') 2 F 2 p + a' 2 )\ 

hi = S22 = /4 (1 + (27ra') 2 HB 2 ) 1 (l - (2W) 2 F 2 p + o' 2 ) 

§33 = 4 ^ (1 + (2na') 2 HB 2 )^ (l - (W) a i£ + ag»)T 

^ = /if P 2 (1 + (2na') 2 HB 2 )^ (l - (2W) 2 F 2 + ^. 

(3.18) 

# a (a; = 1,2,3) are the angular variables on S 3 and Q a p is the metric on the 
unit 3-sphere. The surface of the presumably apparent horizon, which is 
defined locally as a surface whose area variation vanishes along the null ray 
which is normal to the surface, is given by 



J \i=l a=l 



fi 7 H* / g 2 + (27ra') 2 P 6 (l + (27ra') 2 HB 2 ) 

M V (i+^F) 



(3.19) 



where this should satisfy 

dV\ dt= . dz = 0. (3.20) 

In order to simplify the calculations we obtain the equation for the apparent 
horizon using the derivative of the square of (3.19) and get 



2(1 + a ( ' 2 



o ) k 99 - 9 2 (JpM + P) + H ^{f^) 
~ (2W) 2 p 6 (j^M + p)+H^ ^ - 3H^ p~^ (3.21) 
+ (2 7 ra')YHB 2 ( ^(a a' + p) - H^^_ + 3^) = 0, 



cr, 



where g = d x -g. We will use this equation in the following subsections to 
obtain the thermalization time observed on the boundary after we set the 
functional form of g(x~). The calculations have been done for the small and 
general values of the magnetic field, separately. 
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3.1 Weak Magnetic Field Limit 

In this subsection we are going to study the meson thermalization time for 
the small values of constant external magnetic field, i.e. 2ira'B C 1. As it 
was mentioned in section 2 one can set <x (p) approximately to zero for weak 
magnetic field. This will simplify the calculations significantly. The equation 
(3.21) reduces to 

±R*p- 1 gg - Ag 2 + 2(2vra / ) 2 p 6 - 2(2™') 4 p 2 R 4 B 2 = 0, (3.22) 
or equivalently 

X-2n 2 z 4 B 2 2 

-V + 4.-W-0, (3.23) 

where g(x~) = g max y(x~) and g max was introduced in (3.7). In order to 
determine the thermalization time we have to specify the form of the function 
y(x~). Introducing different appropriate functions for y(x~) lead to various 
forms of baryon injection into the system. This function increases form zero 
to a constant g ma x during baryon injection and remains constant after that. 
A suitable function is [12, 14] 

D, x~ < 

g(x~) = { g max ux~, 0<x~<^ (3.24) 

9maxi 7T *^ *^ 

where ux~ is a dimensionless combination. The time-dependent injection of 
baryons leads to the horizon formation on the probe brane where its location 
can be obtained from (3.23). The results are 

3z JX + 16n 2 7i^z s cu 2 - 2n 2 z 4 B 2 1 

t = T ± ~ t^r. • t<z + ^' < 3 - 25a > 



where 



3072n|j7r 4 



u = J B 12 -4608S 6 n|vr 2 A + 24576 1 144n|vr 4 A 2 



+v / 6n|vr 3 A 3/2 x /3456n|j7r 2 A - S 6 ] . 



(3.26) 
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Figure 3: The location of the apparent horizon for t < z + 1/u in the t — z 
plane for {2na')B = 0.09. The red curve represents (3.25a) and the blue line 
shows the earliest tangent null ray. 



In figure 3 the red curve represents the location of the apparent horizon 
for t < z+l/u, (3.25a) 4 . The thermalization time observed on the boundary 
is defined as the time where the earliest null ray tangent to the curve reaches 
the boundary. This has been shown by the point B in figure 3 where the null 
ray (blue line) is tangent to the curve at point A. Hence in the light cone 
coordinates the thermalization time is 

tf h = om = t A + z A . (3.27) 

It is evident from the above relations that the value of the thermalization 
time depends on the external magnetic field. In figure 4(left) the difference 
between the thermalization time in zero and non-zero magnetic field (the 
blue curve), 5t t h = tf h — t^ =0 , is numerically plotted as a function of the 
external magnetic field. It is clear that Stth becomes more and more neg- 
ative as the magnitude of the magnetic field increases. In other words the 
thermalization of mesons or meson dissociation happens more rapidly in the 
presence of the magnetic field. Interestingly the magnetic field dependence 

4 In order to plot this figure we have set the parameters as [12,14]: rig — 34, A = 10, 
w = 10 and R — 1. We will use these numbers throughout the paper to plot different 
figures. 
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of the thermalization time is precisely fitted with 



Hh 



29 x 10" 6 (2W) 2 J B 2 



h B=0 

J th 1 



(3.28) 



where tf h =0 ~ (^2) 1/8 [12,14]. The equation (3.28) reveals that 5t th is 

proportional to a' 2 B 2 (the red curve in figure 4(left)). The decrease in the 
thermalization time is explained by comparing the energy of the system with 
and without the magnetic field. For a static solution the energy is equal to 
the negative of the action evaluated in that solution. Therefore using (2.6) 
the energy density, E = At7 v 3 vo^3) J Jt > ^ or a con fig ura tion with nonzero 
magnetic field is 



and hence 




B= + (3-29) 



SE = E — Eqj 



(3.30) 

where Eq is the energy density when 5 = 0. The equation (3.29) shows 
that the magnetic field increases the value of the energy density. This fact is 
numerically shown in figure 4 (right). Therefore a solution with a larger value 
of the magnetic field has a larger value of energy density and by injecting 
the energy into the system it would be easier for mesons to dissociate. Hence 
t£ < if, 2 for B x > B 2 . 




B 0.030 




Figure 4: Left: Stth for t < z + ^ is plotted in terms of B in the weak magnetic 
field limit. The red curve shows the fitted function. Right: The energy of 
the massless embedding in terms of B in the weak magnetic field limit. 
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Figure 5: Left: The location of the apparent horizon for t > z + 1/u in 
the t — z plane for (27ra')B = 0.09. The black, red and blue lines represent 
(3.25b), t = z + - and the earliest null ray which reaches the boundary. 
Right: 5t t h for t > z + - is plotted in terms of B in the weak magnetic field 
limit. The red curve shows the fitted function. 



Depending on different values of the parameters, the thermalization may 
occur after the baryon injection has been ceased i.e. t > z + 1/w. Then 
the equation describing the location of the apparent horizon on the flavour 
D7-brane is given by 

!6n|^ -V = 0. (3.31) 

The solution to this equation has been previously found in (3.25b). In figure 
5 (left), t = z + l/u and (3.25b) are shown by black and red lines, respectively. 
Similar to what we discussed before the point B represents the earliest time 
that a boundary observer sees the apparent horizon formation on the probe 
brane (point A). The blue line shows an outgoing null ray form the point A 
to the boundary observer. Therefore the thermalization time is 

+ B=0.09 _ + , 

Hh ~tA + ZA, 

1 (3.32) 

tA = Z A + -, 
UJ 

where the value of za is given by (3.25b). The same analysis done before to 
see how the thermalization time changes with magnetic field can be applied 
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here. In figure 5 (right) the thermalization time has been plotted with respect 
to the magnetic field and can be fitted by 

fg = (l - 28 x 1(T 6 (2 W) 2 5 2 ) 4=°, (3.33) 

1 /6 

where tff~° = (^t) [12,14]. Interestingly we observe that although tff~° is 
different in (3.33) and (3.28) the dependence of the thermalization time on 
B is the same for both cases. 




Figure 6: dtth for t < z + - is plotted in terms of B for general values of B. 
The red curve shows the fitted function. 



3.2 General Magnetic Field 

In the previous section we investigated weak magnetic field limit in which 
we could assume ctq(p) ~ 0. Now we would like to generalize it to arbitrary 
values of B where we have to use (3.21) and try to solve it numerically. We 
follow exactly the same steps as the previous section and the result for the 
thermalization time with respect to the value of the magnetic field has been 
shown in figure 6. The blue dots show our numerical results. It can be fitted 
by (red curve in figure 6) 

tf h =(l-4x 10- 5 (2 7 ra') 2 5 L435 )tfr°. (3.34) 

Comparing this with (3.28) and (3.33) we observe that as we increase the 
magnitude of the magnetic field the power of B in the fitted equation deviates 
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from the power 2 of B in those equations. Note that similar to the weak 
magnetic field limit the thermalization time decreases as one increases the 
magnetic field. 

Concluding Remarks 

• The main result is that the meson thermalization(dissociation) happens 
more rapidly in the presence of external magnetic field. A reasonable 
argument for this observation is to evaluate the energy of the system. 
One can see from figure 4 (right) that the energy of the system increases 
as we increase the magnetic field on the probe brane. The figure in- 
cludes only small values of the magnetic field but this argument can 
be generalized to more general values of it. Therefore we can conclude 
that as we inject the energy to the system by the baryon injection the 
mesons can dissociate more easily for higher values of the magnetic 
field. 

• In weak magnetic field limit Stth oc B 2 and for general case Stth oc 
£1.435 Therefore we can see the inclusion of larger values of the external 
magnetic field on the probe brane will lower the power of B in the 
relation for Stth- 

• We have observed that the thermalization time decreases with the raise 
in the magnetic field. This is true up to a specific value of the magnetic 
field which seems to be around 20 and then this behaviour reverses. It 
seems that very large values of the magnetic field distorts the system in 
a way that meson dissociation has already happened due to the large 
energy that the magnetic field induces to the system. Similar observa- 
tion has been done in superconductors where there is a critical value of 
the magnetic field after which the superconducting phase doesn't hap- 
pen. Also the phase transition to the superconducting phase happens 
faster as the magnetic field is raised [20]. 
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